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Aperture synthesis techniques are increasingly being employed to provide high angular resolution images in 
situations where the object of interest is in the near field of the interferometric array. Previous work has showed 
that an aperture synthesis array can be refocused on an object in the near field of an array, provided that the 
object is smaller than the effective Fresnel zone size corresponding to the array-object range. We show here 
that, under paraxial conditions, standard interferometric techniques can be used to image objects which are 
substantially larger than this limit. We also note that interferometric self-calibration and phase-closure image 
reconstruction techniques can be used to achieve near-field refocussing without requiring accurate object range 
information. We use our results to show that the field of view for high-resolution aperture synthesis imaging of 
geosynchronous satellites from the ground can be considerably larger than the largest satellites in Earth orbit. 

OCIS codes: 100.3175 Interferometric imaging; 280.4991 Passive remote sensing; 

http://dx.doi.org/10.1364/XX.99.099999 


The techniques of aperture synthesis were initially de¬ 
veloped for astronomical imaging, where the distance 
to the object being studied could in most cases be as¬ 
sumed to be infinite. This “far field” approximation is 
not valid for more recent applications such as microwave 
remote sensing [T] , terahertz imaging [5] and tracking of 
space debris|3], because the objects are sufficiently close 
to the interferometric array that the curvature of the 
wavefronts from the target at the location of the imag¬ 
ing array cannot be neglected. 

This “near held” condition occurs when the distance R 
between the array and the object is given by i? < 
where Smax is the maximum interferometric baseline 
used and A is the observing wavelength. To give a prac¬ 
tical example, obtaining 10 cm-resolution images from 
the ground of satellites in geosynchronous orbit requires 
an interferometer with a maximum baseline length of 
-Bmax ~ 350 m, assuming an operating wavelength of 
A = 1 fj,m. For such an interferometer, near-held ef¬ 
fects become important for objects closer than about 
10® km, so all geosynchronous satellites, which orbit ap¬ 
proximately 36,000 km from the Earth’s surface, are well 
inside the near held. 

Carter [3] showed that the wavefront curvature oc¬ 
curring in the near held could be compensated for in an 
interferometer by a procedure analogous to refocussing a 
telescope. However the analysis assumed that the max¬ 
imum extent X^ax of the object satishes X^ax ^ 
where R is the distance from the array to the object, 
equivalent to the object being smaller than the hrst 
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Fresnel zone radius for propagation over a distance R. 
This restriction is quite limiting in many cases. In the 
geosynchronous satellite example given above, only ob¬ 
jects much less than 6 m in size satisfy Carter’s assump¬ 
tion. Many satellites in geosynchronous orbit are signif¬ 
icantly larger than this, with “bus” sizes of 10-15 m and 
solar panels up to 50 m in span. 

Lazio [5] discusses a method for getting around this 
restriction using an adaptation of the “w-projection” al¬ 
gorithm used for far-held but wide-angle interferometric 
measurements. This procedure requires the use of ad¬ 
ditional data processing steps and cannot directly use 
standard astronomical image reconstruction software. 

In this letter, estimates are derived for the held of view 
which can be imaged using standard interferometric im¬ 
age reconstruction software with little or no modihca- 
tion. It is shown that objects signihcantly larger than 
the held of view suggested by Carter can be straight¬ 
forwardly imaged in a number of situations of practical 
interest. 

The geometric model for the imaging system is shown 
in Figure [2 Both the interferometer and the object un¬ 
der observation consist of a set of points in arbitrary 
three-dimensional arrangements. In the case of the ob¬ 
ject these points represent the locations of emitters of ra¬ 
diation and in the interferometer these points represent 
collectors of radiation (telescopes at optical wavelengths 
and antennae at radio wavelengths). It is assumed that 
the collectors are clustered around an array center O and 
that the emitters are clustered around a “phase center” 
P. The vector from O to P is denoted by the vector 
R and its length R is called the “range”. The choice of 
the locations O and P is to some extent arbitrary but 
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where '4’iBj) is the instantaneous complex wave ampli¬ 
tude measured at collector j and angle brackets denote 
averaging over a time much longer than the coherence 
time of the radiation. In an optical interferometer these 
visibilities would be obtained by formation of interfer¬ 
ence fringes and measurement of the fringe parameters, 
whereas in a radio interferometer the direct product of 
the measured field amplitudes would be formed and av¬ 
eraged. 

Considering the radiation from a single emitter i of 
strength A^, then using a scalar wave approximation, 
the complex wave amplitude received at collector j is 
given by 


^ A,exp[ 27 ri|f,,|/A] 

m ) « -pg— 


( 2 ) 


where is the vector from collector j to emitter i. This 
vector is given by 


T j^j — 13 j -t- I? -t- 


( 3 ) 


Fig. 1. Geometric model of the interferometer and object. 


a good choice for each location is one which minimises 
the distance to the furthestmost point in the respective 
cluster. 

The vector from P to emitter i is given by Xi = 
{xi,yi, Zi)R where the coordinate system is a Cartesian 
one such that the z axis is parallel to R. The vector from 
array point j to O is given by Bj = (uj,Vj,Wj)R using 
the same coordinate axes. Note that the u, v and w 
coordinates defined in this way differ from the standard 
aperture synthesis variables of the same name because 
they are normalised by the range R rather than by the 
wavelength A. 

An important assumption is that the maximum di¬ 
mension of the object Almax and the maximum dimen¬ 
sion of the interferometer i?max are much less than R. 
In other words |xip -|- -I- <C I and -I- \vj\'^ + 

\wj\'^ 1 for any i or j. This “paraxial” condition is 

more restrictive on the collector array size than the ge¬ 
ometry adopted by Carter [1] who assumed that the col¬ 
lectors can be distributed over a half-sphere surrounding 
the emitters. However, paraxial conditions are satisfied 
in many conditions of practical interest, with the geosyn¬ 
chronous satellite example serving as one case in point. 

The emitters are assumed to be isotropic emit¬ 
ters and incoherent with one another. Only quasi- 
monochromatic radiation at a single wavelength A is 
considered. The interferometer is modelled as an ab¬ 
stract device which measures complex “fringe visibil¬ 
ities” (more accurately “coherent fluxes”, since these 
have units of flux, rather than being normalised to unity 
for a point source) given by 

v,u = {i>mr{Bu)), ( 1 ) 


and so writing R = (0,0, l)i?, the distance rg- = jr^j 
between the emitter and collector is given by 

Tij = {\xi \yi -\- -|- [1 -1- Zi -|- ^ R- 

( 4 ) 

Under paraxial assumptions the denominator in equa¬ 
tion ([^ can be approximated as « R. Substituting 
equation ([^ into equation Q gives an expression for the 
complex fringe visibility 


Vjk oc 


I A, p exp 


2TTid{Xi,Bj,Bk)/X 


( 5 ) 


where d{Xi,Bj, Bk) is the optical path difference (OPD) 
given by 

d(^Xi^ Bj^ Bk) = Tij Tik- ( 6 ) 


Without loss of generality, we can consider the visibil¬ 
ity on a baseline where one collector is situated at 
O and define the OPD as a function of two coordinates 

diX,,B,) = d{X,,Bj,0). (7) 


The following results can be straightforwardly extended 
to any collector pair in the array, since 

diX,, Bj,Bk) = d{X,, B,) - d(X„ Bk). (8) 


Defining a scale parameter a which is the greater of 
Xjnax/R and Bjnax/R, the OPD can be expanded to sec¬ 
ond order in a to give 

d{Xi,Bj) = R (^XiUj + yiVj + Wj + -I- ^Vj + 0(a^)^ 

( 9 ) 

The first two terms in the expansion can be recognized as 
giving rise to the linear dependence of the fringe phase 
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on the off-axis angle of the emitter, identical to that 
for far-held interferometry. The following three terms 
are independent of the emitter position. The w term is 
normally removed using an internal delay in the inter¬ 
ferometer which “points” the array at the phase center. 
The quadratic terms arise from the spherical nature of 
the wavefront from the emitter and so are negligible in 
the far held when 


i?(u| -b v]) < A. (10) 

This is equivalent to the standard far-held condition 

( 11 ) 

In the near-held, the observed visibility phase can be 
corrected to the value which would be observed if the 
source were in the far held by subtracting a phase oh- 
set which is quadratically dependent on the distance be¬ 
tween the collector and O. This correction is equivalent 
(to second order) to the “focussing” correction suggested 
by Carter [3]. We choose to write this phase correction 
in terms of the visibility phase which would be observed 
for a point emitter at the phase centre, so that the phase- 
corrected visibility is given by 


Vjo' = Vjo exp —27rid(0, Bj)/\ 


( 12 ) 


For a collection of incoherent sources then the phase- 
corrected visibility will be given by 


oc ^ X! 1^*1^ 27rid'(Xi, Bj)/X 

where d' is the corrected OPD given by 

d'(Xi, B,) = d(X„ Bj) - d(0,4). 


(13) 


(14) 


If a continuous emitter brightness distribution I(x,y,z) 
is approximated by a set of discrete emitter brightnesses 
given by \Ai\^ oc I(xi,yi, Zi) dxdydz then we recover 
the usual two-dimensional Fourier relationship between 
the phase-corrected visibility and the (depth-integrated) 
brightness distribution 


measured visibility is corrupted by arbitrary “antenna- 
dependent” phase errors (for example due to atmo¬ 
spheric seeing) and solve for these errors simultaneuosly 
with solving for the image intensity distribution. It is 
straightforward to demonstrate that the near-held focus 
correction as written in equation 111 is in an antenna- 
dependent form and so the closure phase/self-calibration 
process will solve for the sum of the near-held phase er¬ 
ror and any phase errors due to seeing. Thus if closure 
phase or self-calibration imaging software is used, then 
near-held visibility data can be used directly without 
requiring pre-processing to provide focus correction. 

Using self-calibration or closure phase imaging soft¬ 
ware to correct the visibility phases has the additional 
advantage that a precise estimate of the range to the 
target is not required in order to derive a precise focus 
correction. For geosynchronous satellite imaging, the 
satellite range must be known to better than 0.01% in 
order to compute an a-priori phase correction accurate 
to 1 radian at 1 micron wavelength, and this kind of 
range accuracy may be difficult to achieve. 

In practice, the high-precision phase correction af¬ 
forded by self-calibration or closure phase software may 
need to be accompanied by lower-precision OPD ad¬ 
justments in hardware in order to deal with temporal 
coherence effects. For example, the fringes formed in 
an interferometer operating at a central wavelength of 
1 /im with a bandpass of 10 nm will have a fringe en¬ 
velope which extends over approximately ±100 ^m of 
OPD. When imaging a satellite in geosynchronous or¬ 
bit using baseline lengths of 350 m, the OPD error due 
to near-held effects will be about 1.7 mm, and so the in¬ 
strumental OPD will need to be adjusted compared with 
the far-held case in order to see high-contrast fringes. 
This hardware adjustment needs only to be accurate at 
the sub-100 /rm level: any remaining phase errors will 
be antenna-dependent and so can be taken care of by 
self-calibration or closure phase software. 


The imaging equation given in equation (15) is correct 
to second order in the scale parameter a. 
to third order in a gives 


Ixpanding d' 


d'(Xi, Bj) = R [xiUj + yiVj - e + 0{a'^)] (16) 

where e is a third-order error term given by 


V'qOC J exp[2ni{UjX+ Vjy)] I{x,y,z)dz^ dxdy, 

(15) 

where Uj = Ruj/X and Vj = Rvj/X are the standard in¬ 
terferometric “u-v” coordinates. In other words, provid¬ 
ing the visibility data is phase-corrected for focus, then 
images can be reconstructed from near-held interfero¬ 
metric data using standard far-held aperture synthesis 
software. 

It is important to note that explicit phase correction 
of the visibility data is not required if closure phase 
or equivalent “self-calibration” [5] image reconstruction 
methods are used. These methods assume that the 


e = {x,Uj ± y,Vj ) {zi ± ) ± ^ Wj {xl ± y ■) ± ^ (u^ ±). 

(17) 

It is notable that all the terms in e contain a factor which 
depends either on the depth of the array Wj or the depth 
of the object Zi. This means that if both the array and 
object are planar so that wj = Zi = 0 for all i and j then 
the third-order terms vanish and only fourth-order error 
terms need to be considered. 

In general however, the array and/or the object will 
be three-dimensional, and in this case the field-of-view 
limitations will arise when the third-order terms lead 
to phase errors which are comparable to a radian, i.e. 
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when Re ~ A. The terms in e are all products of object- 
dependent distances and array-dependent distances. If 
the interferometer is large compared to the object, the 
largest terms in e will be of order , while 

if the object is larger than the interferometer, then the 
largest terms will be of order Bmax^max/^^- 

In the former case, the third-order phase errors will 
be negligible when 




< R 


R 


(18) 


The factor R g^ is approximately the size of a reso¬ 
lution element Xres in tbe interferometric image, so the 
maximum size of object which can be imaged is approx¬ 
imately i?/i3max resolution elements across. Thus the 


field of view expression in equation (181 can also be writ¬ 
ten 




<C Xr 



(19) 


maximum baselines needed to provide the desired res¬ 
olution, and so instead equation (21) must be used to 
estimate the field of view. This gives X^ax ~ 1.9 km 
which is still much larger than the 50 m dimensions of 
the largest satellite solar panels. 

Geosynchronous satellite imaging at microwave fre¬ 
quencies may be less straightforward: equation (19) sug¬ 
gests that, if 10 cm resolution is required at a wavelength 
of 1 cm, third-order phase aberrations can become be¬ 
come significant for fields of view of as small as 1 meter. 
However, the field of view improves quadratically as the 
resolution requirement is relaxed: if 50 cm resolution is 
adequate for a given investigation, then imaging a field 
of view of 25 m would be unproblematic. 

The above analysis shows that standard self¬ 
calibration imaging algorithms should work well for 
geosynchronous satellite imaging applications at visi¬ 
ble and near-infrared wavelengths. In this case, factors 
other than the imaging algorithm, such as providing suf¬ 
ficient dense sampling of the u — v plane [7] , are likely to 
be the key challenges to successful imaging of the largest 
objects. 


If instead the object is larger than the interferometer, 
then the phase errors will be negligible when 

( 20 ) 

\ -^max / 

which can be written in terms of the interferometric res¬ 
olution Xj-es as 


^max \/^res7?. 


( 21 ) 


Equations ( |I^ and give expressions for 

the field of view for imaging using standard self¬ 
calibration/closure phase image reconstruction methods. 
It is interesting to apply these results to some cases of 
practical interest. 


A one-dimensional teraherz aperture synthesis exper¬ 
iment is described by Bandyopadhyay et al. [5] using an 
array of collectors with Hmax = 5 cm and R = 40.9 cm. 
The spatial resolution of this arrangement at the op¬ 
erating frequency of 0.535 THz (corresponding to A = 
0.542 mm) is x^es ~ 4.43 mm, so using equation (19) 
gives a value for the field of view of about 8 resolution 
elements or 36 mm. Bandyopadhyay et al. image a point 
source approximately 10 resolution elements away from 
the optical axis and find the image reconstruction ade¬ 
quate, indicating that the field of view given by equa¬ 
tion (19) may be pessimistic in this case. 

In the case of a geosynchronous satellite observation 
aiming for a;res=10cm resolution at A = 1 /rm, then the 
maximum object size that can be observed is given by 
equation (19) as 10 km. This is much larger than the 
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